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1.0  INTRODUCTION 

1.1  Jlbjective 

(C)  "^To  develop  a combined  track-classify  algorithm  based  on  state-of-the-art  target 
motion  analysis  and  space-time  classification  methods.  To  map  the  performance  envelopes  of 
the  algorithm  relative  to  the  AN-SQQ-23  sonar  system  parameters  and  other  active  systems 
as  time  and  manpower  permit.  To  write  a proposed  technical  approach  (PTA)  to  active  classi- 
fication in  contemporary  sonar  systems  based  on  these  performance  envelopes. 


1.2  Background 

(C)  A recent  report  by  D.  G.  Olson  (Reference  1)  outlines  rough  performance  pre- 
dictions for  STARLITE  classification  method  based  on  detection  performance  predictions  of 
the  AN-SQQ-23  sonar  system  and  the  Weikhorst  geometric  constraint  equations  for  the  STAR- 
LITE classification  method.  However,  the  Weikhorst  constraints  are  incomplete  because  they 
do  not  take  echo-to-noise  ratio  into  account.  It  is  a fairly  safe  assumption  that  an  operational 
STARLITE  processor  would  be  sensitive  to  echo-to-noise  ratio  which,  in  turn,  is  dependent 
on  sonar  system  parameters,  and  track  geometry.  The  task  which  remains  is  to  design  a 
specific  operational  procedure  for  using  STARLITE  information  in  combination  with  other 
sonar  outputs  in  an  optimal  fashion,  and  to  measure  rather  than  predict  its  performance. 


1.3  Approach 

(C)  A specific  classification  procedure  using  both  STARLITE  and  target  motion  in- 
formation is  to  be  designed  and  tested.  It  will  be  designed  using  the  AN-SQQ-23  PAIR  sonar 
system  parameters  as  a basis  for  defining  available  sensor  information.  It  will  be  tested  using 
computer  simulation  techniques,  incorporating  both  computer  generated  target  echoes  and 
actual  sonar  echo  sequences  from  the  AN-SQQ-23  PAIR  system.  The  artificial  echo  data  is 
necessary  to  maintain  parameter  control  for  measuring  performance  over  a wide  range  of  pa- 
rameter values.  Tape  recorded  sonar  echoes,  although  limited  to  a fairly  restricted  range  of 
parameter  values,  is  necessary  to  lend  the  necessary  degree  of  credibility  to  the  results. 

(C)  Most  simply  stated,  target  track  information  will  be  used  to  “predict”  STARLITE 
processor  outputs.  The  predicted  outputs  will  be  compared  with  measured  STARLITE  outputs 
to  determine  whether  or  not  a sufficient  amount  of  self-consistency  exists  to  justify  making  a 
contact  classification.  The  rationale  for  combining  target  motion  information  with  STARLITE 
clues  follows  simply  from  the  fact  that  the  STARLITE  processor  output  is  heavily  dependent 
on  relative  track  geometry.  If  the  track  geometry  can  be  determined  through  independent 
measurements  (range,  bearing,  Doppler),  much  of  the  ambiguity  associated  with  a STARLITE 
processor  output  can  be  removed.  One  of  the  basic  problems  to  be  resolved  is  to  determine 
how  many  echo  returns  are  necessary  to  gain  a sufficiently  accurate  track  solution  based  on 
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the  range,  bearing,  and  Doppler  measurement  accuracy  of  the  sonar  system  in  question.  Once 
this  basic  comparison  scheme  is  perfected,  other  classification  clues  may  be  added  to  the  de- 
cision rule  to  help  improve  performance. 

1.4  Plans,  Milestones 

(C)  Figure  1 is  a program  flow  chart  in  which  specific  tasks  are  indicated  by  rectangles 
and  milestones  by  octagonals.  A milestone  is  defined  as  a point  in  the  program  flow  at  which 
the  decision  must  be  made  to  continue,  to  go  back  and  modify,  or  to  stop  and  redefine  the 
program.  Normally,  these  are  shown  following  the  completion  of  each  task.  This  does  not 
mean,  however,  that  the  task  must  be  carried  through  to  the  end  before  such  a decision  can  be 
reached.  An  NUC  Technical  Note  will  be  written  at  each  milestone  point  to  provide  a tangible 
basis  for  program  evaluation. 

(C)  At  the  beginning,  the  flow  diagram  consists  of  four  main  branches.  These  are: 

I.  Data  production  and  collection  for  performance  mapping. 

II.  Development  and  performance  mapping  of  a STARLITE  processor  for  the 
AN-SQQ-23  sonar  system. 

III.  Target  motion  algorithm  selection,  modification,  and  performance  mapping. 

IV.  Display  development  for  data  monitoring  and  operator  clue  evaluation. 

These  four  branches  converge  to  a single  branch  in  which  an  algorithm  is  assembled  to  use 
both  STARLITE  and  TMA  information  to  classify  the  contacts.  The  output  of  this  algorithm 
is  to  be  both  an  automatic  classification  decision  and  a graphical  data  display  for  operator 
clue  evaluation. 

The  following  text  is  a list  of  milestones  and  explanations. 

(C)  Milestone  IA  - Computer  Generated  Echoes.  At  this  point,  all  software  programs 
for  supplying  computer  generated  echoes  with  tape  recorded  reverberation  will  be  complete 
and  debugged.  Effects  due  to  variable  target  range,  speed,  aspect,  turning  rate,  own  ships  mo- 
tion and  beam-forming  will  be  incorporated  into  these  echoes.  The  tape  recorded  reverbera- 
tion will  be  taken  from  USS  Glennon  dual  dome  array  (PADLOC  II).  This  data  will  be  filtered 
through  digital  filter  in  the  ASDAC  System. 

(C)  Milestone  IB  - Tape  Recorded  Echoes.  At  this  point  all  tape  recorded  echo  data 
from  the  AN-SQQ-23  PAIR  receiver  will  be  in-house  and  evaluated  for  quality;  i.e..  whether  or 
not  spurious  frequencies  or  other  effects  such  as  overclipping  arc  present  that  would  seriously 
impair  data  usage. 

(C)  Milestone  HA  - STARLITE  Software  Completion.  At  this  point  software  will  be 
complete  and  debugged  for  a STARLITE  processor  algorithm.  Preliminary  performance  checks 
will  also  be  completed  to  insure  that  the  programs  are  efficiently  composed.  Since  this  algor- 
ithm will  later  be  used  for  many  repeated  runs  in  a monte-carlo  type  performance  evaluation, 
it  is  essential  that  it  use  as  little  computer  time  as  possible. 
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(C)  Milestone  IIB  - STARLITE  Processor  Performance  Mapping.  At  this  point  a map- 
ping of  performance  envelopes  of  the  STARLITE  processor  will  be  completed  for  both  a wide 
range  ot  parameter  values  with  computer  generated  echoes  and  restricted  range  of  parameter 
values  with  tape  recorded  echoes  from  the  AN-SQQ-23  PAIR  sonar  receiver.  A tentative  pa- 
rameter list  includes  target  range,  bearing,  aspect,  speed,  turning  rate,  own  ships  motion,  and 
echo-to-noise  ratio.  It  would  be  impossible  to  investigate  all  combinations  of  parameter  values 
due  to  computer  time  limitations.  Therefore,  for  the  computer  generated  input  data,  only 
certain  areas  of  parameter  space  consistent  with  expected  performance  levels  of  the  PAIR  re- 
ceiver for  specific  track  geometries  will  be  investigated.  These  geometries  will  be  designed, 
based  on  the  most  common  types  of  target  encounters,  as  discussed  in  Reference  2.  The  pa- 
rameter space  for  the  tape  recorded  input  data  is,  of  course,  determined  by  the  parameter 
values  available  in  the  tape  library. 

(C)  Milestone  III  A - TMA  Survey  and  Preliminary  Performance  Evaluation.  At  this 
point  a series  of  monte-carlo  type  computer  simulations  will  be  completed  which  will  demon- 
strate the  target  motion  analysis  capabilities  of  a Kalman  tracking  algorithm  with  inputs  avail- 
able from  the  AN-SQQ-23  sonar  system.  Specifically,  the  ability  to  estimate  target  course  and 
speed  as  a function  of  ping  number  will  be  investigated  for  seven  simple  track  geometries  taken 
from  Reference  1.  The  results  of  this  study  will  be  used  to  identify  critical  parameters  in  the 
course  estimation  problem.  Once  the  critical  parameters  are  identified,  the  decision  can  be 
made  whether  to  modify  this  algorithm  to  handle  them  more  effectively  or  to  seek  another 
candidate  algorithm  for  this  purpose.  A tentative  goal  is  to  be  able  to  estimate  target  course 
within  ±10°  in  5-  10  pings  for  target  less  than  10  Kyds  in  range. 

(C)  Milestone  IIIB  - TMA  Performance  Mapping.  At  this  point  performance  mapping 
ot  the  selected  TMA  algorithm  will  be  completed  for  the  same  parameter  set  as  chosen  for  the 
STARLITE  processor  performance  mapping. 

(U)  Milestone  IVA  - MSDD  Display  Interface.  At  this  point  the  MSDD  display  on 
loan  to  ASDAC'S  will  be  interfaced  directly  with  the  CDC  1 700  computer  in  ASDACS,  bypas- 
sing the  UNIVAC  667  computer.  The  steps  necessary  to  reach  this  point  arc:  to  build  line 
drivers  and  a cable  to  bypass  the  ASDACS  digital  patch  panel  for  the  output  channel  from  the 
C DC  1700  to  the  display;  to  adapt  and  check  out  the  existing  digital  input  channel  through 
the  delogger  interface  in  ASDACS;  to  write  display  driver  software  routines  for  CDC  1700. 

Once  these  steps  are  complete,  desired  display  format  will  be  easily  programmable. 

(C)  Milestone  IVB  - Display  Formats.  At  this  point,  display  formats  will  be  completed 
and  debugged  which  will  display  STARLITE  processor  outputs  and  predicted  outputs  generated 
by  the  I MA  algorithm.  These  formats  will  be  used  for  comparison  of  human  operator  classifi- 
cation for  the  selected  classification  vector  against  an  automatic  classification  performed  by 
the  STARL1TE-TMA  algorithm  in  the  next  task  area. 

(C)  Milestone  V - Program  Convergence.  At  this  point,  1 B,  IIB.  IIIB,  and  IVB  must 
be  completed  so  that  the  combined  STARLITE  - TMA  classification  algorithm  can  be  assem- 
bled and  evaluated  as  a single  entity. 
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(C)  MILESTONE  VI  - Track -Classify  Algorithm  Evaluation.  At  this  point  the  com- 
bined STARL1TE  - TMA  classification  algorithm  will  be  evaluated  on  an  ROC  curve  basis, 
averaged  over  the  track  geometries  used  previously  for  STARLITE  and  TMA  performance 
maps.  The  ROC  analysis  will  be  completed  for  both  operator  and  automatic  due  evaluation. 
A tentative  design  goal  is  to  classify  at  90%  with  a false  alarm  rate  of  10%  within  three  pings 
after  acquiring  a target  track. 


1.5  Current  Status  Summary 

(C)  A simulated  sonar  echo  generation  algorithm  has  been  completed  for  a stationary 
target  and  constant  velocity-receiving  platform.  Equations  for  a nonstationary  target  and  re- 
ceiving platform  with  yaw,  surge,  and  sway  velocity  components  have  been  derived. 

(C)  The  time  and  frequency-domain  implementations  of  a STARLITE  processor 
have  been  compared.  The  frequency-domain  processor  appears  superior  when  target  bearing 
is  inaccurately  known  and  when  matched  filtering  is  desirable  for  processing  gain.  The  time- 
domain  processor  has  the  advantage  of  being  easily  implemented  and  performs  well  if  target 
bearing  could  be  accurately  determined. 

(C)  Tape-recorded  reverberation  has  been  obtained  for  use  with  synthetically  gen- 
erated echoes  and  has  been  quantitatively  evaluated. 

(U)  A digital  display  has  been  partially  interfaced  directly  with  the  CDC  1700  com- 
puter in  ASDACS.  F urther  software  refinements  are  needed  for  general  usage. 


2.0  AREAS  OF  EFFORT 
2.1  Applied  Research 

2.1.1  Echo  Simulation 

The  objective  ot  this  effort  is  to  develop  a computer  program  to  generate  simulated 
submarine  echoes  for  a dual-receiving  array  active  sonar  system.  This  echo  generation  algor- 
ithm will  sum  the  echoes  from  individual  reflectors  on  the  target  model  (Reference  3),  as 
well  as  duplicate  this  spatial  beamforming  of  the  SQQ-23  PAIR  beamformers  (Reference  4). 
The  simulated  echoes  will  include  the  effects  of  the  following  own  ship  and  target  motion 
components,  own  ship’s  roll,  yaw,  surge,  sway,  and  constant  velocity  along  own  ship’s 
heading:  and  target  constant  velocity  and  constant  turning  rate.  The  above  motion  compo- 
nents will  be  resolved  into  equivalent  relative  motion  between  each  hydrophone-reflector 
pair  to  compute  the  resultant  time  delay  and  Doppler  time-scale  factor  for  each  echo  com- 
ponent. 
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(C)  The  target  model  reported  in  Reference  5 has  been  modified  to  include  aspect- 
dependent  shading,  documented  in  Reference  3.  The  introduction  of  shading  vastly  im- 
proved the  target  model,  as  is  evidenced  by  Figures  2 4.  Figure  2 shows  the  target  strength 
pattern  of  the  unshaded  model,  and  Figure  3 shows  that  of  the  shaded  model.  Figure  4 
shows  a measured  target  strength  pattern  for  the  USS  Nautilus  (Reference  6).  Note  that  the 
difference  between  the  target  strength  at  beam  aspect  and  that  at  bow  or  stern  aspect  is  10  dll 
tor  the  Nautilus,  9 dB  lor  the  shaded  model,  and  6 dB  for  the  unshaded  model.  These  com- 
parisons  verify  not  only  that  shading  vastly  improves  the  target  model,  but  also  that  the  shaded 
model  is  a very  reasonable  two-dimensional  acoustic  representation  of  a submarine. 

(C)  In  addition  to  the  target  model,  the  derivation  of  the  defining  equations  for 
the  echo  simulation  algorithm  has  been  completed  and  documented  in  Reference  7.  The 
equations  define  exactly  the  echoes  received  at  both  the  forward  and  aft  arrays  of  a dual 
array  sonar  system  for  (a)  the  case  of  a stationary  target  and  quasi-stationary  sonar  plat- 
form and  (b)  the  case  of  a nonstationary  target  and  sonar  platform.  The  equations  for  the 
quasi-stationary  case  give  the  echoes  as  functions  of  the  parameters:  range,  bearing,  aspect 
angle,  own  ship  constant  velocity,  and  the  orientation  of  the  beamformer.  Those  for  the 
nonstationary  case  give  the  echoes  as  functions  of  range,  bearing,  aspect  angle,  own  ship 
constant  velocity,  orientation  of  the  beamformer,  own  ship  constant  velocity,  target  con- 
stant turning  rate  and  constant  velocity,  displacement  due  to  roll,  period  of  roll,  maximum 
angular  displacement  due  to  roll,  displacement  due  to  yaw,  period  of  yaw,  maximum  angu- 
lar displacement  due  to  yaw,  sway  velocity  and  displacement,  and  surge  velocity  and  dis- 
placement. 

(U)  The  remaining  task  in  the  evolution  of  the  echo  simulation  algorithm  is  a sen- 
sitivity analysis  of  the  defining  equations  of  the  simulated  echoes  to  determine  which  short- 
cuts can  be  taken  in  the  computation  of  the  simulated  echoes,  and  thereby  minimize  the 
complexity  of  the  algorithm  and  its  software  implementation. 

2.1.2  Starlite  Processors 

(C)  The  ultimate  goal  of  this  effort  is  to  determine  the  optimal  STARLITE  re- 
ceiver lor  the  AN/SQQ-23  sonar  system  relative  to  performance  envelopes  and  processor 
complexity.  The  initial  effort  has  been  to  compare  time-  and  frequency-domain  processors 
(Figure  5)  relative  to  four  criteria: 

1.  Mean  absolute  error  in  aspect  angle 

2.  Standard  deviation  of  aspect  angle  error 

3.  Mean  output  correlation  peak  width  at  0.9  maximum  peak  height 

4.  Average  main-to-secondary  peak  height  ratio 

5.  Average  main  peak  height. 

The  results  are  summarized  in  Table  1. 
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(A)  Time-domain  implementation  of  STAR  LITE  (assumes  the  transmit  signal  is  a linear  FM  slide). 


OPTIONAL 


(B)  Frequency  domain  implementation  of  STARLITE  (assumes  the  transmit  signal  is  'band-limited'). 


(C)  Figure  5.  Time-  and  Frequency-domain  Implementations  of 
STARLITE  Processors. 
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(C)  TABLE  1.  STATISTICAL  COMPARISON  OF 
TIME-  AND  FREQUENCY-DOMAIN  PROCESSORS 


Time  Domain 

Frequency  Domain 

. . - - - - - - 1 

Mean  |error|  of  a 

0.12° 

0.856° 

/\ 

Standard  deviation  a 

1.06° 

3.0° 

Mean  main  correlation  peak  height 

0.647 

0.919 

Mean  main  correlation  peak  width 
(0.9  down.  0.707  no  work) 

| 

5.84  Hz 
(equivalent) 

9.975  Hz 

Mean  main-to-secondary 
correlation  peak  ratio 

1.88 

1.215 

(C)  These  statistics  were  obtained  by  computing  the  output  of  the  processors  for  an 
ensemble  of  36  computer  generated  echo  plus  noise  records.  The  echoes  were  generated  from 
the  stationary  target  model  discussed  in  sections  2.1.1.  The  parameter  values  for  these  trials 
were:  R = 5 Kyd,  aspect  angle  = 45°,  and  target  bearing  of  90°.  The  sonar  array  geometry  and 
transmit  pulse  waveform  parameters  were  those  of  the  PAIR  sonar  receiver  using  a 161  milli- 
second FM  slide  at  5 Khz  center  frequency.  Examples  of  the  outputs  of  the  time-  and 
frequency-domain  processors  are  shown  in  Figure  6a  and  6b. 

(C)  These  results  were  computed  assuming  the  target  bearing  was  known  exactly. 
However,  the  time-domain  approach  is  much  more  sensitive  to  errors  in  bearing  measurements 
than  is  the  frequency-domain  processor.  For  example,  for  a bow  or  stern  target  at  a range  of 
10  Kyds  and  relative  bearing  between  30°  and  1 50°,  the  aspect  angle  error  can  be  as  high  as 
27°  for  a bearing  error  of  3°.  For  the  frequency-domain  processor,  the  error  is  less  than  2° 
for  the  same  parameter  values.  Therefore,  the  frequency-domain  processor  is  preferable  unless 
the  bearing  error  could  be  held  to  less  than  0.5°. 


2.1.3  Tracking  Study 

A small  study  contract  has  been  awarded  to  Singer- Librascope  Co.  to  compute  expected 
error  variances  in  target  course  estimates  as  a function  of  ping  number  after  detection  for  the 
Singer- Librascope  tracking  algorithm.  Estimates  of  course  from  a track  solution  can  be  used  to 
resolve  ambiguity  in  STARL1TE  processor  outputs.  The  objective  of  this  study  is  to  determine 
how  many  echoes  are  needed  to  obtain  a sufficiently  accurate  track  solution  given  the  predic- 
ted accuracies  of  the  active  range  and  bearing  estimates  of  the  PAIR  receiver.  The  test  tracks 
selected  for  analysis  were  selected  from  Reference  1 in  order  to  be  consistent  with  previous 
efforts  within  Code  606  to  analyze  STARL1TE  processor  performance.  It  is  expected  that 
once  a contact  is  classified,  the  STARL1TE  outputs  could  be  used  to  improve  a fire  control 
solution. 
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(C)  f igure  6(b).  Simulated  Output  C'orrelogram  for  the  Frequency- Domain 
Implementation  of  STARLITE:  <p  ~ 45°,  R - 5 Kyd.  = 40°,  S-N-R  = 0 dB. 
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2. 1  4 Detector  Performance  Comparison 

(C'l  This  work  was  performed  in  response  to  a special  request  by  SHIPS  (OOV1C)  that 
certain  correlation  detector  performance  curves  in  Reference  X be  recomputed  for  a false 
alarm  rate  of  .01.  The  performance  curves  presented  here  are  for  RDT  short  pulse  CW  trans- 
missions from  the  AN/SQS-23  sonar  receiver  documented  in  Reference  8.  The  received  echoes 
were  from  a beam  aspect  target.  The  FM  transmission  data  also  documented  in  Reference  8 
was  not  available  for  this  effort.  The  noise  weighted  linear  correlator  (NWLC).  simple  linear 
correlator  (SLC).  and  polarity  coincidence  correlator  (PCC)  were  compared  relative  to: 

1 . Their  theoretical  performance  curves  ( Figure  7 ). 

2.  Their  performance  in  an  ambient  sea  noise  limited  environment  (Figure  8). 

3.  Their  performance  in  a reverberation  limited  environment  (Figure  9). 

(U)  In  Figure  7 the  curve  for  the  NWLC  was  computed  assuming  an  ideal  stationary 
bandlimited  noise  background.  The  curve  for  the  SLC  was  computed  assuming  an  ideal  white 
noise  spectrum.  The  curves  for  the  PCC  were  also  computed  for  a white  noise  spectrum  but 
with  false  alarm  rates  of  .03 1 and  .006  instead  of  .01 . Due  to  the  limited  output  dynamic 
range  of  the  PCC.  the  values  .03 1 and  .006  were  the  closest  usable  values  bracketing  01  The 
curve  for  the  PCC  for  a false  alarm  rate  of  .01  would  lie  approximately  half  way  between  those 
for  .031  and  .006. 

(U)  The  conclusions  suggested  by  Figure  7 are  the  same  as  derived  in  Reference  8 for 
the  higher  false  alarm  rate  of  0.1.  The  best  obtainable  performance  in  stationary  bandlimited 
noise  is  roughly  4 dB  below  the  best  obtainable  performance  in  white  noise,  and  approximately 
4 dB  is  lost  by  using  the  PCC  instead  of  the  SLC  in  white  noise. 

(U)  Figure  8 shows  the  measured  performance  curves  for  the  NWLC,  SLC.  and  PCC  in 
ambient  sea  noise  for  a false  alarm  rate  of  .01 . Here,  again,  bracketing  curves  were  computed 
tor  the  PC  ( . Since  the  narrowband  filtered  ambient  sea  noise  is  very  similar  to  stationary  band- 
limited  noise,  the  NWLC  performs  best  since  it  is  designed  for  this  type  of  noise  background. 
Figure  9 shows  the  measured  performance  curves  for  a reverberation  limited  environment. 
Again,  curves  for  obtainable  false  alarm  rates  bracketing  .01  were  computed  for  the  PCC'. 

Here  it  should  be  noticed  that  although  reverberation  is  bandlimited.  the  SLC  outperforms  the 
NWLC  by  3 dB  Furthermore,  the  PCC  is  the  best  receiver  for  S/N  less  than  I dB.  These  results 
are  apparently  due  to  the  non-stationary  nature  of  a reverberation  limited  noise  background. 


2.2  Research  Support 

2.2.1  Data  Collections 


2. 2. 1.1  Current  Status  — (C)  The  data  collection  plan  outlined  in  the  previous 
progress  report  consisted  of  two  phases  which  have  subsequently  been  greatly  reduced  in 
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(O  Figure  7.  Theoretical  Performance  Cunes  for  the  Noise  Weighted  Correlator  (Ci  I. 
Simple  Linear  Correlator  (C2 ) and  Polarity  Coincidence  Correlator  (C3)  for  a PFAofO.Ol 
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(C)  Figure  8.  Empirical  Performance  Curves  for  the  NWLC  (C, ).  SLC  (C2).  and 
PCC  (Cj)  in  an  Ambient  Sea  Noise  Limited  Background  for  a PFA  of  0.01 . 
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scope  because  of  cutbacks  in  funding.  The  objective  of  Phase  I was  to  obtain  analog  beain- 
tornied  reverberation  data  from  the  AN-SQQ-23  PAIR  receiver  A secondary  purpose  was  to 
check  the  quality  of  the  data  before  proceeding  to  the  much  more  expensive  Phase  II 

<(')  In  the  revised  program,  the  reverberation  data  will  not  be  from  the  PAIR  system 
but  rather  from  the  PADLOC  II  system  aboard  the  USS  Glennon.  The  PADLOC  II  arrays  are 
very  similar  to  those  in  the  PAIR.  However,  the  FM  pulse  durations  and  bandwidths  are  slightly 
diflerent,  hence,  the  Glennon  reverberation  data  will  not  be  “matched”  to  the  PAIR  pulse 
parameters.  I his  data  is  now  in-house.  However,  being  twice  dubbed,  it  suffers  from  several 
spurious  effects  outlined  in  section  2.2. 1.2. 

(t ) The  original  objective  of  Phase  11  was  to  obtain  analog  beamformed  echo  data. 
Passive  data,  and  a simultaneously  recorded  digital  tape  containing  own  ships  motion  com- 
pensation information,  sonar  status  words,  and  a time  code  for  matehing  the  echo  data  with 
track  plots  taken  on  the  Saint  Croix  tracking  range. 

(U)  In  the  revised  plan.  Phase  II  will  not  include  the  digital  tapes.  This  will  result  in 
a substantial  savings  in  funds  since  the  necessary  hardware  and  software  for  syncronizing  the 
analog  and  digital  tapes  during  the  dubbing  operation  will  no  longer  be  needed.  However,  it 
will  also  result  in  a substantial  loss  in  desired  information. 

2-2.1 .2  PADLOC  II  Reverberation  Data  Evaluation  (C)  Representative  portions  of 
5 reels  of  analog  beamformed  reverberation  data  from  the  PADLOC  II  array  aboard  the  USS 
tilennon  were  aurally,  visually,  and  analytically  evaluated  for  use  in  future  development  work. 

I hese  tapes  are  dubs  of  dubs  and  contain  several  spurious  effects  which  must  be  avoided  in 
future  use  of  these  tapes. 

(C)  The  analog  beamformed  data  from  channels  recorded  on  both  the  fore  and  mid 
domes  on  the  PADLOC  II  array  were  played  back  through  a high  fidelity  speaker  system  for 
aural  detection  of  undesirable  signals.  This  data  was  also  digitized  at  20  kHz  and  plotted  for 
visual  inspection.  An  analytical  evaluation  was  performed  by  computing  an  FFT  over  a data 
window  of  0. 1 second  once  every  hall  second  for  each  ping  cycle.  The  resulting  sequence  of 
l l-fs  was  plotted  lor  each  ping  cycle  to  locate  spurious  frequency  component  throughout 
the  ping  cycle.  See  Figures  10  and  1 1. 

(C)  Visual  inspection  shows: 

1.  All  returns  are  clipped  for  two  to  lour  seconds  after  the  ping  cycle  begins.  As 
would  be  expected,  longer  transmit  signals  result  in  longer  periods  of  clipping. 

I his  creates  a dead  area  in  target  range  in  which  analysis  of  data  cannot  be 
performed. 

2.  A few  sequences  are  clipped  over  the  entire  ping  cycle  making  them  useless. 

3.  Adjacent  beams  of  the  loreward  array  are  attenuated  by  different  amounts. 

I his  may  have  occurred  during  one  of  the  tape  dubbing  procedures. 
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(C)  Figure  10.  Upper  Graph:  Amplitude  Versus  Time  Plot  of  Digitized  Samples  from  Beam  9 of  Fore  Dome 
of  PADLOC  II  Array  for  One  Ping  Cycle  of  Reverberation  for  a 240  msec  400  Hz  FM  Upslide  Centered  at 
5 kHz.  Lower  Graph:  Time  Sequence  of  Fast  Fourier  Transforms  for  a 0. 1 sec  Data  Window  Taken  Every 
0.5  Seconds. 
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(C)  Figure  1 1.  Upper  Graph:  Amplitude  Versus  Time  Plot  of  Digitized  Samples  from  Beam  8,  Mid  Dome 
of  PADLOC  II  Array  for  One  Ping  Cycle  of  Reverberation  for  a 1 20  msec  400  Hz  FM  Upslide  Centered  at 
5 kHz.  Lower  Graph:  Time  Sequence  of  Fast  Fourier  Transforms  for  a 0.1  sec  Data  Window  Taken  tverv 
0.5  Seconds. 
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4.  Bottom  reverberation  appears  to  be  present  in  many  of  the  returns  8 to  1 0 
seconds  after  commencement  of  the  ping  cycle  on  most  of  the  data.  The  re- 
verberation appears  somewhat  more  “echo  like”  than  the  volume  reverberation. 

(C)  Frequency  analysis  indicates: 

1.  Background  noise  levels  are  45  to  54  dB  down  from  the  largest  frequency  com- 
ponent in  the  passband.  This  is  satisfactory. 

2.  Spurious  frequency  component  spikes  are  present  at  0.6  kHz.  1 kHz.  6.5  kHz. 
and  7 kHz.  These  spikes  are  generally  1 2 to  24  dB  above  the  background  noise 
level.  These  spikes  can  be  eliminated  by  proper  filtering  of  the  data. 


2.2.2  ASDACS 

2.2.2. 1 Current  Status 

2.2. 2. 1.1  General  Purpose  Software  (U)  The  following  routines  are  general  purpose 
FORTRAN  compatible  programs  and  subprograms  written  for  the  ( DC  1 700  computer. 

1.  (U)  PROGRAM  STATIS:  Computes  the  mean,  variance,  probability  density 
function  of  a given  ensemble  according  to  a specified  resolution.  The  density  function  and 
the  distribution  function  are  given  in  tabular  form  and  are  optionally  plotted.  (The  program 
uses:  FUNCTION  XMEAN.  FUNCTION  VAR.  SUBROUTINE  INPUT,  SUBROUTINE  BIN 
and  various  Calcomp  plotting  routines.) 

2.  (U)  SUBROUTINE  HOROUT  (X.Y.N)  and  SUBROUTINE  RD1N1  (X.Y.N.IS. 
IU).  These  are  general  purpose  utility  routines  that  facilitate  handling  of  floating  point  arrays 
between  different  computers.  (Used  here  to  transfer  information  between  CDC  1604  and 
CDC  1 700  computers  and  these  routines  are  used  on  both  machines.)  HOROUT  is  used  to 
write  an  array  of  N values  of  X vs.  Y on  a particular  tape  unit.  The  first  record  of  a section 
of  data  contains  the  value  of  N as  well  as  the  minimum  and  maximum  values  of  the  array  Y. 
Subsequent  records  contain  eight  values  of  X followed  by  eight  values  of  Y alternately  until 
N values  of  each  are  written.  The  tape  format  is  in  FORTRAN  E-notation  suitable  for  listing. 
RDIN1  is  used  to  read  in  the  data  in  the  same  format,  i.e.,  N values  of  X and  Y from  unit  IU 
and  data  section  specified  by  IS. 

3.  (U)  SUBROUTINES  CON VOUT  (1X.IY.N),  CONV,  PACK.  OUTBUF.  These  are 
a set  of  utility  routines  written  for  the  CDC  1604  computer  to  produce  binary  tapes  that  are 
directly  readable  by  the  CDC  1 700.  Included  is  a conversion  routine  to  convert  CDC  1604 
real  numbers  to  CDC  1 700  formatted  real  numbers. 

4.  ( U)  SUBROU I INE  EN VGEN  ( F.N.FH.F L.I  S.E ).  This  subroutine  generates 
the  envelope  E of  N real  values  ol  data  F given  FH  (the  high  frequency  component).  FL  (the 
low  frequency  component)  and  FS  (the  sampling  frequency).  (This  routine  is  also  operational 
on  the  CDC  1604  computer.) 
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5.  (U)  FUNCTIONS  XMAX1F  (X,N),  XMINIF  ( X,N).  ABMAXF  (X,N):  generate 
the  maximum,  minimum,  and  absolute  maximum  of  N values  of  real  array  X respectively. 

6.  (U)  FUNCTION  FRF1  1 (X.N):  generates  the  error  function  integral  of  X and 
it  indicates  N,  the  number  of  terms  used  in  the  exponential  expansion.  The  error  function  is 

e 1 dr  and  given  X this  function  generates  F accurate  to  about  10 


7.  (U)  FUNCTION  ERFINV  (E):  computes  the  inverse  error  function,  i.e.,  given 
E this  function  returns  X (see  the  integral  above)  accurate  to  about  10  3 . 

8.  (U)  PROGRAM  BC'DPLT:  a general  purpose  plotting  routine  that  accepts  BCD 
formatted  tape  (via  RD1N1 ) of  real  numbers,  scales  the  data  by  a given  scale  factor  and  op- 
tionally plots  several  sections  of  data  on  one  set  of  axes. 

9.  (U)  SUBROUTINE  1MPLT  < X,Y.N,XST,X1NC,XLN.YST.Y1NC.YLN ).  This 
subroutine  plots  an  impulse  plot  of  N points  of  X vs.  Y on  the  Calcomp  plotter.  The  plot  = 
XLN  by  YLN  inches  originating  at  point  ( XST,  YST)  with  XINC  and  YINC  as  the  X and  Y 
increments  per  inch,  respectively. 

10.  (U)  SUBROUTINE  C1RPLT  (X.Y.N.XYMX).  This  subroutine  produces  a polar 
plot  of  10  X 10  inches  with  maximum  Y value  being  XYMX.  It  assumes  that  the  X array  is  an 
angle  measurement  given  in  radians.  The  routine  performs  the  transformation  from  cartesian 
to  polar  coordinates  and  produces  the  plot.  For  the  use  by  this  routine,  the  standard  Calcomp 
AXIS  subroutine  was  rewritten  in  order  to  allow  axis  divisions  to  be  stepped  off  from  maxi- 
mum to  minimum  over  the  length  of  the  axis. 

1 1.  (U)  PROGRAM  AMPSEC:  computes  the  Fresnel  integral.  This  program  was  used 
to  check  the  Fast  Fourier  Transform  of  a waveform. 

12.  (U)  SUBROUTINE  CORREL  ( F.G.NF.NG.R ):  correlates  real  array  F with  real 
array  G and  stores  the  result  in  array  R.  NF  is  the  number  of  values  of  F.  NG  is  the  number  of 
values  of  G,  and  it  is  assumed  that  R is  dimensioned  to  NF  + NG  — 1 . 

13.  (U)  FUNCTION  CORRELF  (TAU.G.TG.F.TF.NG.NF).  This  function  assumes 
to  have  arrays  of  NG  values  of  TG  vs.  G and  NF  values  of  TF  vs.  F and  computes  one  step  of 
a correlation  of  F with  G assuming  a displacement  of  amount  TAU  of  TG  relative  to  TF.  The 
function  value  is  the  value  of  the  correlation.  (This  function  was  also  implemented  on  the 
CDC  1604  computer.) 

14.  (U)  SUBROUTINE  NEWTON  (AC'.N.R.RO):  computes  the  root  R of  a given 
function  F using  Newton’s  method  with  RO  as  the  initial  estimate  and  to  within  specified 
accuracy  AC.  N is  the  number  of  iterations  it  takes  to  achieve  the  result.  It  is  assumed  that 
FUNCTIONS  F and  FPR  are  given.  F as  the  function  and  FPR  its  derivative. 
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15.  ( U)  SUB  KOI'  I INE  SOK  I ( X.k.N ):  sorts  an  array  X of  N real  numbers  in  as- 
cending order,  k is  a key  integer  index  array  that  indicates  the  relative  position  of  elements 
in  the  unsorted  array  so  that  arrays  parallel  to  X can  be  sorted  by  simply  addressing  elements 
indirectly  through  k. 

16.  (U)  FUNCTIONS  IRDB  ( lO.k.lU)  and  IWTB  ( lO.k.lU).  These  are  FORTRAN 
callable  read  and  write  routines,  respectively.  They  read/write  k binary  words  of  array  IO 

1 rom/onto  unit  1U  where  1U  is  a magnetic  tape  unit.  The  function  value  is  an  indication  of 
tape  unit  status  at  the  completion  of  the  operation. 

17.  (U)  SUBROUTINFs  BREAD  (1U.10.L.1FL)  and  BWRITE  (lU.lO.L.lFLl 
These  routines  are  the  generalization  of  IRDB  and  IWTB  respectively  and  are  designed  to 
include  input/output  from/onto  the  disk  as  well  as  the  tape  units.  L is  a three-word  parame- 
ter array  and  IFL  is  the  status  indicator. 

18.  (U)  FUNCTION  MOTION  (IU,ICODE):  allows  FORTRAN  programs  to  have 
complete  motion  control  of  magnetic  tape  units  such  as  rewind,  backspace,  advance  to  end- 
of-file.  etc.  The  function  value  is  an  indication  of  tape  unit  status  at  the  end  of  the  operation. 

19.  ( U)  PROGRAM  NOISE:  generates  a specified  number  of  points  of  white  noise 
(zero  mean,  unit  variance)  on  the  CDC  1604  computer  with  the  aid  of  the  standard  random 
number  generator  and  writes  them  on  magnetic  tape  in  standard  FORTRAN  E-format  so  that 
the  white  noise  is  accessible  to  the  CDC  1 700  computer. 

20.  (U)  PROGRAM  NARBNOIS  was  written  for  the  CDC  1604  computer  to  gen- 
erate narrow  band  noise  by  using  white  noise  and  cascading  it  through  a triple  stage  digital 
filter.  The  format  of  output  is  the  same  as  for  the  white  noise  (see  above). 

21.  (U)  SUBROUTINE  MATINV  (A. N,B.M,DET).  This  is  an  adaptation  of  a sub- 
routine to  the  CDC  1 700  machine  that  was  written  for  the  CDC  1604  computer.  The  subrou- 
tine performs  the  inversion  ol  the  N X N matrix  A with  accompanying  solution  of  M columns 
of  linear  equations  (optional).  DET  is  the  determinant  of  the  matrix  A. 

2. 2. 2. 1.2  Digital  Display 

A.  HARDWARE: 

1 (U)  A 32  circuit  0. +6  volt  to  0.  1 2 volt  output  level  shifter  and  cable 

driver  has  been  built  to  adapt  the  CDC  1 700  digital  output  voltage  levels  to  the  voltage  levels 
required  by  the  Modified  Sonar  Digital  Display  (MSDD)  data  and  control  lines. 

2.  (U)  An  input  level  shifter,  originally  intended  to  adapt  time  and  sonar 
parameter  logging  code  translator  data  to  the  voltage  levels  required  by  the  ASDACS  digital 
inputs,  already  exists  and  will  be  used  for  input  data  from  the  MSDD. 

3.  (U)  Work  remaining  to  be  done  on  the  ASDACS-MSDD  INTERFACE 

includes: 

(a)  Building  Input  and  Output  cables  to  run  between  ASDACS  and  the  MSDD; 
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(b)  Including  a self-contained  power  supply  in  the  output  level  shifter  and 
cable  driver; 

(c)  Acquisition  and  installation  of  logic  boards  for  two  additional  digital  output 
words  to  be  committed  solely  to  the  MSDD. 

B SOFTWARE: 

1.  (U)  The  following  subroutines  have  been  written  and  together  make  up  a 
Fortran  callable  display  generator  and  driver. 

(a)  SETBUF  (1BUF.N):  initializes  a display  buffer.  1BUF.  N words  long  in  the 
calling  routine. 

(b)  DLINE  (X.Y.DX.DY ):  draws  a vector  starting  at  point  (X.Y ) with  X deflec- 
tion DX.  Y deflection  DY. 

(c)  CLINE  (DX.DY):  using  the  end  point  of  the  last  vector  as  a starting  point, 
draws  a vector  with  X deflection  DX,  Y deflection  DY. 

(d)  EVESET  (XADJ.IEV.ISIZ);  sets  certain  parameters  used  by  the  display  to 
determine  symbol  size  and  type,  and  modify  X position. 

(e)  STEXT  (X.Y.ICODE.IB);  draws  one  of  42  alphanumeric  or  special  symbols 
determined  by  ICODE,  at  one  of  three  brightness  levels  determined  by  IB, 
at  a position  (X.Y).  Symbol  size  (full  or  half  size  available)  is  set  by  parame- 
ter ISIZ  in  EVESET. 

(f)  SPADE  (X.Y.IAMP.1DOP):  draws  a sonar  event  symbol  at  positions  (X.Y) 
where  symbol  type,  length,  and/or  brightness  is  determined  by  IAMP  and 
IDOP,  modified  by  parameter  IEV  in  EVESET. 

(g)  P1CTUR:  performs  the  cyclic  output  of  the  display  buffer  generated  by  the 
above  subroutines  to  the  MSDD. 

2.  (U)  Test  routines  to  checkout  the  above  software  and  the  MSDD  and  inter- 
face hardware  have  also  been  written. 

3.  (U)  General  purpose  software  work  remaining  to  be  done  includes: 

(a)  A routine  to  perform  input  of  data  from  the  MSDD; 

(b)  Extension  of  the  routine  STEXT  to  provide  the  capability  of  outputing  a 
message  based  on  a Fortran  format  statement  or  an  integer  or  floating  point 
number; 

(c)  An  executive  routine  to  display  a function  list  and.  on  the  basis  of  inputs 
from  the  MSDD.  select  one  of  several  processing  and  display  routines. 
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2.3  Auxiliary  Tasks 

2.3.1  Travel 

(U)  Messrs.  G.  Turton  and  J.  Teeter  visited  ARL.  Sperry  Corp..  and  SHIP  OOVIC  to 
discuss  procedures  for  obtaining  tape  dubs  of  AN/SQQ-23  PAIR  beamformed  data  during  the 
week  of  18  November.  1969. 


2.4  Technical  Reports  Issued 

1.  Turton.  G.A..  "Problem  Summary  for  SF-1 1-121-106/8132  for  the  Period  from  1 Jan- 
uary 1969  to  30  June  1969  (U),”  NUC  TN  283.  September.  1969  (CONFIDENTIAL). 

2.  Teeter.  J.L..  "A  Two  Dimensional  Target  Model  <U),"  NUC  TN  315.  10  December 
1969  ifttNI miiVII  \l  i 
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